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The axial- vector Ward-Takahashi identity is used to derive mass formulae for neutral pseudoscalar 
mesons. Flavour symmetry breaking entails non-ideal flavour content for these states. Adding that 
the r\ is not a Goldstone mode, exact chiral-limit relations are developed from the identity. They 
connect the dressed-quark propagator to the topological susceptibility. It is confirmed that in the 
chiral limit the jj' mass is proportional to the matrix element which connects this state to the vacuum 
via the topological susceptibility. The implications of the mass formulae are illustrated using an 
elementary dynamical model, which includes an Ansatz for that part of the Bethe-Salpeter kernel 
related to the non-Abelian anomaly. In addition to the current-quark masses, the model involves 
two parameters, one of which is a mass-scale. It is employed in an analysis of pseudoscalar- and 
vector-meson bound-states. While the effects of SU(Nf — 2) and SU(Nf = 3) flavour symmetry 
breaking are emphasised, the five-flavour spectra are described. Despite its simplicity, the model is 
elucidative and phenomenologically efficacious; e.g., it predicts rj-n' mixing angles of ~ —15° and 
7r°-?7 angles of ~ 1°. 
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I. INTRODUCTION 

Flavour symmetry breaking has long been of interest. 
For example, it showed up in the application of current 
algebra to strong interaction phenomena. In QCD no 
two current-quark masses are equal. Isospin (SU (2)- 
flavour) breaking is determined by the current-mass dif- 
ference m„ — rrid, while S'?7(3)-fiavour breaking can be 
measured via m s — (m u + m £ ;)/2. The c- and 6-quark 
current-masses are too large for any sensible discussion 
of larger flavour symmetry groups but the light-quark 
mass differences also have an impact on the spectrum of 
hadrons containing a heavy-quark. 

It is of interest to explore and determine the effect of 
these differences in current-quark mass throughout the 
hadron spectrum. This leads one to consider the dif- 
ference in mass between charged and neutral hadrons. 
Part of that splitting is electromagnetic in origin but 
constraining the strong component is necessary before 
one can know just how large that electromagnetic contri- 
bution might be. 

We focus herein on the strong interaction component 
alone. To be specific, we concentrate on exploring the ef- 
fect of flavour-symmetry breaking on pseudoscalar- and 
vector-meson masses. These sectors are of particular in- 
terest because any reliable calculation of pseudoscalar 
meson masses must involve a consideration of the axial- 
vector Ward-Takahashi identity. Moreover, since when 
viewed simply vector mesons are spin-flip partners of 
the pseudoscalars, it is natural to examine how the l - - 
0~ mass splitting evolves with current-quark mass and 
mass difference. On the other hand and in addition, one 



might find that the strong breaking effects in mesons can 
be interpreted judiciously and used to inform results for 
baryons. 

In our analysis we employ the Dyson-Schwinger equa- 
tions (DSEs), pedagogical introductions to which can be 
found in Refs. [ij, [2j. The approach is particularly well 
suited to the study of bound states upon which sym- 
metries and the manner in which they are broken have 
a heavy impact [3|. Within this framework, using the 
rainbow-ladder truncation - the lowest-order in a sys- 
tematic symmetry-preserving truncation scheme 0, Hj], 
effects of flavour-symmetry breaking were reported in 
Ref. [a] . Those results can provide a useful compari- 
son with ours. However, while in one sense a simpler 
interaction is employed herein, we have a different per- 
spective and will demonstrate effects that arise in pro- 
ceeding beyond the leading-order truncation. This is 
especially noteworthy in connection with neutral pseu- 
doscalar mesons, for which the non-Abelian anomaly 
plays an important role p}. 

In Sec. [II] we present the U(Nf) axial-vector Ward- 
Takahashi identity in its general form. It necessarily 
includes a contribution from the non-Abelian anomaly, 
which is explicated. We elucidate a couple of exact re- 
sults that follow from the spectral feature that the rf 
mass is much larger than that of other light-quark pseu- 
doscalar mesons. In Sec. lIIII we derive and discuss exact 
mass formulae for pseudoscalar mesons. Section IIVI in- 
troduces a model that enables the illustration of impli- 
cations of these formulae. It also reports a calculation 
of the masses of ground state pseudoscalar and vector 
mesons for Nf — 5, and covers the phenomena of mix- 
ing between the 7r°, r\ and rf . The results enable us to 
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provide estimates of the non-electromagnetic component 
of the neutron-proton mass difference and the masses of 
hitherto unseen B*j mesons. Section [V] recapitulates on 
the main qualitative features emphasised by our study. 



II. AXIAL- VECTOR WARD-TAKAHASHI 
IDENTITY 

The axial-vector Ward-Takahashi identity is basic to 
any study of pseudoscalar mesons. The impact of this 
statement of chiral symmetry and the pattern by which 
it is broken is felt even by heavy-light |8[ and heavy-heavy 
bound-states @. 

The general form of the identity can be expressed 

P^(k;P) = S-^k+^^ + i^S- 1 ^) 

-2iM ab T\(k-P)~ A a (k-P), (1) 

where P — P\— P2 is the total and k the relative momen- 
tum between the amputated quark legs Eq. ([1]) is 
fully renormalised and it is important that the product 
Ai ab T\ does not mix with other operators under renor- 
malisation. 

In Eq. (P), {F a \a = 0,...,Nj - 1} are the gen- 
erators of U(Nf) in the fundamental representation, 
orthonormalised according to X,xT a T h = \8 ah . The 
dressed-quark propagator S — d\&,g[S u , Sd, S s , S c , S&, . . .] 
is matrix-valued with nonzero entries that can be ex- 
pressed in various equivalent forms, e.g., 



S(k) 



1 



Z(k 2 



i-y ■ kA(k 2 ) + B(k 2 ) i~/-k + M(k 2 )' 



(2) 



The propagator is determined by the gap equation 



SfcO-^Zafo.p + M^+Efr), (3) 
S(p) = ZiJ g 2 D^(p - q)^ lfl S(q)T a v (q,p), (4) 

wherein: represents a Poincare invariant regu- 

larisation of the integral, with A the regularisation 
mass-scale H; is the dressed-gluon propagator; 

r v (q,p) is the dressed-quark-gluon vertex; and M hm — 
diag[m^ m , m^ m , m|? m , to]?™ 1 , nvj! m , . . .] is the matrix of A- 
dependent current-quark bare masses. The quark-gluon- 
vertex and quark wave function renormalisation con- 
stants, Zi^(( 2 , A 2 ), depend on the gauge parameter, the 
renormalisation point, £, and the regularisation mass- 
scale. The gap equation is completed by the renormali- 
sation condition 



s(pY 



ip 2 =C : 



= i-fp + M{C), 



(5) 



where -M(C) is the matrix of renormalised (running) 
current-quark masses whose nonzero entries obey 



with Z4 the Lagrangian mass renormalisation constant. 
In Eq. (fTJ) we have defined 



M ab =trj? [{T a ,M}F b } , 



(7) 



where the trace is over flavour indices. 

The inhomogeneous axial- vector vertex in Eq. |T]) sat- 
isfies 

[n^p)] tu = z 2 [ 151 ,T\ u 

+ f [S(q + )T^(q;P)S(q_)} sr KZ°(q,k;P), (8) 

Jq 

where P — p\ —pi = q\—qi and r,. . . , u represent colour, 
Dirac and flavour indices. The pseudoscalar vertex Tg 
satisfies an analogous equation driven by the inhomo- 
geneity Z^ b T a . 

The final term in the last line of Eq. (fT]) expresses the 
non-Abelian axial anomaly. It can be written 

A a (k;P) =S-\k + )5 a0 A u {k;P)S- 1 (k-), (9) 

with 

Au(k; P) = Jd A xd A y e ^ + —k-y) Nf (pO q ^ g (0) ^ . 

(10) 

Here the matrix element represents an operator expecta- 
tion value in full QCD; the operation in Eq. @ ampu- 
tates the external quark lines; and 

Q{x) = i-^-tvc [e^upaF^Fpaix)] = df.K^x) (11) 

47T 

is the topological charge density operator, where the 
trace is over colour indices and = \\ a F^ v is the 
matrix-valued gluon field strength tensor. It is plain and 
important that only A a=0 is nonzero. NB. While Q(x) 
is gauge invariant, the associated Chern-Simons current, 
K^, is not. 

If one imagines there are Nf massless quarks, then 
dynamical chiral symmetry breaking (DCSB) is a neces- 
sary and sufficient condition for the a / components 
of Eq. ([T]) to guarantee the existence of N 2 — 1 massless 
bound-states of a dressed-quark and -antiquark [3j . 

However, owing to Eq. @, a = in Eq. (fT]) requires 
special consideration. One case is easily covered; viz., 
it is clear that if .4° = 0, then the a = component 
of Eq. (fTJ) is no different to the others and there is an 
additional massless bound-state in the chiral limit. 

On the other hand, the large disparity between the 
mass of the ry'-meson and the octet pseudoscalars sug- 
gests that A ^ in real-world QCD. Let's consider this 
possibility and proceed by allowing that T®^ might pos- 
sess a longitudinal massless bound-state pole. In this 
case one can write 



^r BS (k;P) 



Z 2 (C 2 ,A z )m^(A) = Z i (C 2 ,A A )m(0- 



(6) 



+^°75 P) + 7 • kk ■ PG%(k; P) 

+a^P v H° R (k;P) + f^(k;P) 



(12) 
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where F° , G% H° R and T° 5fl (k; P) are regular as P 2 -> 0, 
P M f^ ~0(P 2 ), r ss (fc;P) is the possible bound-state's 
canonically normalised Bethe-Salpeter amplitude an d r * 
is its residue. The amplitude takes the general form [12j 

T BS (k; P) = 2^° 75 [iE BS (k; P) + 7 ■ PP BS (fc; P) 
+7 • kk ■ PG BS (k; P) + a^kpPuHBsik; P)] . (13) 

Since in these circumstances one can write 

A°(k; P) = _P° 75 [i£ A (k; P) + 7 • PT A {k- P) 

+7 • fcfc • P£U(fc; P) + a^P„H A (k; P)] , (14) 

then the Goldberger-Treiman relations of Ref. [3| become 

2r^P ss (fc;0) = 2P (fc 2 )-^(fc;0),(15) 

F°(k;0) + 2r A F BS (k;0) = A Q (k 2 ) - T A (k; 0), (16) 

G^fc;0) + 2r^G ss (fc;0) = 2^,(fc 2 ) - ^(fc; 0),(17) 

P^(fc;0) + 2^P B 5(fc;0) = -H A {k;0), (18) 

where Aq, Bq characterise the gap equation's chiral limit 
solution. NB. A massless pole in A°(k; P) is incompatible 
with Eq. ([I]). 

It now plain that if 



e A (k;0) = 2B {k 2 ), 



(19) 



then r A E B s(k;0) = 0. This being true, then the 
homogeneous Bethe-Salpeter equation (BSE) also pro- 
duces r° A F BS (k;0) = = r A G BS (k;0) = 2r A H BS (k;0). 
Hence, Eq. (fl9|) guarantees that cannot posses a 
massless pole. The converse is also true; namely, the 
absence of such a pole requires Eq. (|19p . It is noteworthy 
that in the neighbourhood of P 2 = 0, Eqs. (TB) - (18} 
thus provide pointwise relations between A°(k;P), the 
dressed-quark propagator and the regular part T^^k; P). 

Equation (|19[) is a necessary and sufficient condition for 
the absence of a massless bound-state pole in Fg . We 
are discussing the chiral limit, in which case E>o(k 2 ) 7^ 
if, and only if, chiral symmetry is dynamically broken. 
Hence, the absence of an additional massless bound-state 
is only assured through the existence of an intimate con- 
nection between DCSB and an expectation value of the 
topological charge density. 

This noteworthy connection is further highlighted by 
the following result, obtained through a few straightfor- 
ward manipulations of Eqs. (I}, (9} and (fTQ|) : 

,A 

(qq)° c = - lim Z 4 (( 2 , A 2 ) tr CD / S°(q,Q (20) 

^ A^°° Jq 



L / d 4 x(q(x)t l5q ( x )Q(0)) a , (21) 



where here the superscript "0" indicates that the quan- 
tity is calculated in the chiral limit. The absence of a 
Goldstone boson in the a — channel is only guaranteed 
if this explicit identity between the chiral-limit vacuum 
quark condensate and the vacuum polarisation generated 
by the topological charge density is satisfied. 



III. MASS FORMULAE 

A wide range of additional observations are possible, 
some of which are canvassed in Ref. [7j . Herein we will 
derive those that are especially relevant in the context of 
this work. 

Equation (T} is an identity that connects two and three 
point functions in QCD. It applies at all values of the to- 
tal momentum P, in particular, at the location of bound- 
state poles. To exploit this we extend Eq. (fT2"]) and ob- 
serve that in the neighbourhood of such a pseudoscalar 
pole, whether massless or massive, 



rL(pi,P2) 



fa p 



P 2 +ml 



ir a 5 (pi, P2 )\ p2 



P 2 + TO; 



r,r,(fc;P) 



(22) 



P 2 



. r ffi (fc;P) 

7V;, 

'(k;P); 



(23) 



viz., each vertex in Eq. (fTJ) is expressed as a simple pole 
plus terms regular in the neighbourhood of this pole, with 
T Wi (k; P) representing the bound-state's canonically nor- 
malised Bethe-Salpeter amplitude [Hj], where i = la- 



bels the lightest pseudoscalar bound-state, i 
next lightest, and so on. In Eqs. ([22} and (|2"3"|) 



f a P 

<(C) 



Zotl 



Z 4 tr 



Jq 



75 Xtt, (q; P) 



1, the 

(24) 
(25) 



where Xm(k;P) = S(k+)T Vi (k; P)S(fc-), k± = q±P/2. 
These residues are gauge invariant and cutoff indepen- 
dent. NB. The nature of r A in Eq. (fT2"|) is now clear. 

While there is certainly no bound-state pole in the in- 
verse of the dressed-quark propagator, the opposite can 
be true of the term associated with the topological sus- 
ceptibility; namely, we must consider 



\P 2 +m 2 r .^0 



-T m (Je;P) 



P 2 +m 2 i 
+ A 0rc Hp llP 2), (26) 



where 



(27) 



Using Eqs. (22}, {23]) and (25} in the axial- vector Ward- 
Takahashi identity, we arrive at a mass formula for pseu- 
doscalar mesons: 



mU:=2M ab pl 



(28) 



It is valid for current-quark masses of any magnitude. 

For nondiagonal mesons this is naturally the same for- 
mula as derived in Refs. d, [H. Allowing for the fact 
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that the Standard Model requires observable particles to 
be eigenstates of the electric charge, it yields, e.g., with 

T K+ = J7i_ 

m 2 K+ f K + 
and similarly, 

■m 2 D+ f D + 



[m u (0 +m s (()]p K+ (() 



[m d (0+m c (0}p D +(0- 



(29) 



(30) 



(NB. With our normalisation, f K + = 113 MeV experi- 
mentally.) Again, these formulae are valid for arbitrarily 
large, or small, current-quark masses. The Gell-Mann- 
Oakes-Renner relation is a small quark mass corollary 
[H, El and aspects of their implications for mesons con- 
taining a heavy-quark, or two, are detailed in Refs. [1,0]. 

A novelty of Eq. (I28p is its validity for charge neutral 
mesons. For example, in the case Nf = 3 one derives for 
the neutral pion 



f„0 













r n 3 1 

Prr" 


f 8 
Jit 







+ 


M 3x3 




Ptt" 


fO 

L Jtt° - 












I Ptv° J 



(31) 



where 



3x3 



WHO 



3 "H-10 



3 mi_io 



mn_2 



3 "li_io 
2 Wn_2 



§ m m 



, (32) 



with m Q ^ 7 = a m u + /3 rrid + 7 m s . 

In the isospin symmetric case; i.e., m u — nid, M 3x 3 
exhibits no mixing between T z and J 70 ' 8 . This signals 
that the flavour content of the ir° is described solely by 
F 3 . Therefore Eqs. ^ and (J2SJ) give / 8 = = f° , 
p^a = = p° , and Eq. ([2^]) yields v v o — 0. Hence, in 
this instance the complete content of Eq. (f3T|) is 



(33) 



mlotfo = [m u (0+m d (0]plo(0- 



This is not true, however, for m v 
subsequently illustrate. 

For Nf = 3 one also obtains 



7^ rrid, as we shall 



(34) 



r f~ 




' " 








~ P%~ 











+ 


M 3x3 




P% 




f° 

L Jrr - 




. n, > . 








.9%. 




r f, ' 













' 9%' 




Jn' 







+ 


Af 3 x3 




P%> 




f, 
L J n - 




. Hr l' . 








.9%' 





(35) 



Naturally, on the domain in which an expansion in 
current-quark mass is valid; viz., m(lGeV) < 50 MeV 
flij . E qs. J3TT) . (|34]) and (|35|) reproduce current algebra 
results fl5j j. 

Of importance is a prediction of the manner by which 
the rf is split from the octet pseudoscalars by an amount 



that depends on QCD's topological susceptibility. This 
is most easily illustrated by considering the U(Nf) limit, 
in which all current-quark masses assume the single value 
m(£). In this case the complete content of Eq. (|35p is the 
statement 



m^ = rv+2m(C)/v(C). 



(36) 



Plainly, the r/ is split from the Goldstone modes so 
long as n v i [l6j]. Numerical simulations of lattice- 
regularised QCD have confirmed the relationship repro- 



duced here 



It is argued [20j, [21( that in QCD 

1 



(37) 



and it can be seen to follow from the gap equation, the 
homogeneous BSE and Eqs. ([24). (|25]) that 



Jn' 



'N r . 



P° V '(0- 



One thus obtains 



7f + 2m(C) 

J ri' 



P°n'(Q 

f°> 
J n' 



(38) 



(39) 



The first term vanishes in the limit N c — > oo while the 
second remains finite. Subsequently taking the chiral 
limit, the rf mass approaches zero in the manner charac- 
teristic of all Goldstone modes. (NB. One must take the 
limit N c — > oo before the chiral limit because the proce- 
dures do not commute (22|. ) These results are realised 
in the effective Lagrangian of Ref. [23| in a fashion that 
is consistent with all the constraints of the anomalous 
Ward identity 0]. 



IV. MESON MASSES: EXEMPLIFYING 
EFFECTS OF MIXING 

A. Model defined 

Implications of the exact results presented above can 
be illustrated and further elucidated by way of a simple 
kernel for the gap and BSEs. We write 



K = Kr 



K A , 



(40) 



where Kl is the leading order in the systematic and 
symmetry preserving truncation explained in Refs. 0>[E]; 
namely, a dressed-ladder interaction: 



{K L )%(q,p;P) = 
-G((p- q ) 2 )D^(p- q ) 









7„y 


ts 


. 7 "T 



,(41) 



wherein D i ^° c (k) is the free gauge boson propagator and 
Q{k 2 ) represents an effective coupling. For the latter we 
use the simple model introduced in Ref. [27j 



g{k 2 ) = (2Tifg 2 k 2 5\k) 



(42) 
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e.g. 



FIG. 1: An illustration of the nature of Ka\ viz., the con- 
tribution to the Bethe-Salpeter kernel associated with the 
non-Abelian anomaly. All terms have the "hairpin" structure 
illustrated in the lower panel. No finite sum of such interme- 
diate states is sufficient. Straight lines denote quarks, with /i 
and fa independent, and springs denote gluons. 



with Q a constant that sets the mass-scale. The model is 
ultraviolet-finite and hence one can remove the regulari- 
sation mass-scale to infinity and set the renormalisation 
constants to one. The infrared enhancement exhibited 
by Eq. (|42|) is sufficient to provide for confinement and 
DCSB, as explained, e.g., in Sec. 2.2 of Ref. @. Moreover, 
in practice it has many features in common with a class 
or renormalisation-group- improved effective- interactions; 
and its distinctive momentum-dependence works to ad- 
vantage in reducing integral equations to algebraic equa- 
tions that preserve the character of the original. There is 
a drawback: the simple momentum dependence can lead 
to some model-dependent artefacts, but they are easily 
identified and hence not generally cause for serious con- 
cern. 

One weakness hampers us, however. The model gen- 
erates Bethe-Salpeter amplitudes whose dependence on 
the relative momentum, k, is unrealistic. In an inter- 
nally consistent definition it is described by S 4: (k). Hence 
one cannot obtain values for the independent overlaps 
and p%. in Eq. l|2"8|) and so the mass formulae can- 
not be directly verified. Nevertheless, we will see their 
imprint in the calculated results for meson masses. Ver- 
ification of Eq. (|28p ispossible with the interactions em- 
ployed, e.g., in Refs. [(| [H, [28[ . Indeed, that has already 
been done for channels in which the rainbow-ladder trun- 
cation is a good approximation, such as charged pseu- 
doscalar mesons and neutral heavy-heavy pseudoscalar 
mesons [1, [Hj]. 

In Eq. (|4"0|) , Ka is a novel addition that we use to 
model effects owing to the non-Abelian anomaly. Its in- 
clusion takes us beyond ladder-truncation and is thus an 
expedient which is a dynamical extension of that em- 
ployed in Ref. [29(. It can be argued from Eqs. (fTU|) and 
that an anomaly-related contribution to a meson's 
Bethe-Salpeter kernel cannot contain external quark or 
antiquark lines that are connected to the incoming lines: 
purely gluonic configurations must mediate, as illustrated 
in Fig.[TJ Furthermore, it is straightforward to see that 



no finite sum of gluon exchanges can serve this purpose. 
Indeed, consider any one such single contribution in the 
chiral limit. It will be proportional to the total mo- 
mentum and hence vanish for P = 0, in conflict with 
Eq. (|36)) . This lies behind the need for something like the 
Kogut-Susskind ghost flq | - (NB. The resummed kernels 
explored in Refs. [3(| l3ll . [33 do not resolve such vac- 
uum polarisation diagrams [5j and thus cannot generate 

As in Ref. [331 ] , with these observations in mind we em- 
ploy 

(K A )%{q,p;P) 
= -f ((g ~ P?) {cos 2 9 ( [C7s]r s fc75]t« 

+ sin 2 9 t [«7 • P^] rs [«7 ' Pis]** } > (43) 
£(fc 2 ) = (2^) 4 £<S 4 (fc), (44) 

where £ is a dimensionless coupling strength. In prin- 
ciple, £(fc 2 ) would also depend on the total momentum 
but for simplicity we ignore that herein. In proposing 
Eq. (|4*3"|) we have also used the fact that Eq. (|4"2"|) only 
supports a pseudoscalar Bethe-Salpeter amplitude of the 
form 



1^(^ = 2.7^75 [igT+l-Pg?] 



(45) 



namely, as described above, the interaction requires the 
constituents' relative momentum to vanish. The angle 
8^ controls the relative magnitude of the two possible 
contributions to the kernel. 

The remaining piece of Eq. (|43[) is the flavour matrix 



diag[ 



1 



1 



1 



Mf = M f (s = 0), 



(46) 
(47) 



where Mt(s) is the mass function for a quark of flavour 
/ [see Eq. ©J. Equation (|4"T|) defines a dynamical 
constituent-quark mass. It differs from the Euclidean 
constituent-quark mass (e.g., Ref. [HI), but is easier to 
calculate and is likewise a renormalisation point invariant 
in QCD. This term introduces a nonperturbative mass- 
dependence, which models that arising from the dressed- 
quark lines that complete a "U-turn" in the so-called 
hairpin diagram in Fig.[T] 



B. Parameters fixed 

The model has three parameters in addition to which 
there are Nt current-quark masses. We determine the 
current-quark masses and Q in Eq. (|42[) by applying the 
model to charged-pseudoscalar and vector meson ground- 
states. Since Ka doesn't contribute in these channels, 
this corresponds to a rainbow-ladder treatment of those 
states, which is plausibly accurate to < 10% [HJ for light- 
quark mesons and becomes precise for heavy-heavy sys- 
tems 31 1 . 
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The rainbow gap equation is obtained from Eq. (|3j) 
with 

Ztu(p) = - I (K L ) t r u s (q 7 p;P)S sr (q) . (48) 
Jq 

(NB. It is plain upon insertion that Ka defined in 
Eg. (|43|) does not modify Eq. (|48f .~) The gap equation 
is solved to obtain the matrix dressed-quark propagator, 
which is then used to complete the homogeneous BSE: 

T H (k;P) = f [S(q + )r H (q + ,q-)S(q-)} S r(K L y t s(q,k;P). 

Jq 

(49) 

Equation (|49[) can viewed as defining an eigenvalue 
problem. There can be and is no mixing between charged 
and neutral mesons so the eigenvector in the case of 
charged pseudoscalars can be written 

r ^0P)= E 2F a K[ip a 1+7 -PpZ}, (50) 

d=l, 2,4,5,6, 7,... 

with the index selecting all Nf(Nf — 1) nondiagonal gen- 
erators of SU(N f ). One inserts Eq. flSOj! into Eq. (gSJ and 
evaluates the spinor trace to arrive at an equation with 
the structure 

P=[K H ,]P, (51) 

with p — columnfpl, pj, pf , p 2 ,, . . .]. (The procedure is 
made explicit, e.g., in Ref. [30].) The matrix Kh 5 has 
Nf(Nf — 1) eigenvalues, {A^ 5 }, and eigenvectors, each 
of which depend on the value of P 2 . One has a solution 
of the homogeneous BSE when one of those eigenvalues 
acquires the value "one" and the mass of the associated 
bound-state is the value of P 2 for which this occurs; viz., 

K ff5 ) 2 = {-P 2 |A^ 5 (P 2 ) = 1}. (52) 

At this point the related eigenvector is that meson's 
Bethe-Salpeter amplitude. If n eigenvectors assume the 
value "one" at the same value of P 2 , then one has an 
n-fold degeneracy. 

With the interaction in Eq. (j42|) . the eigenvector asso- 
ciated with the vector mesons of U(Nt) has the general 
form 

N f 

r k (n = J2 2 Ta b ■ fA < + a ^4 p " «s] > ( 53 ) 

a=l 

where {e*(P)|A = — 1, 0, 1} is the polarisation four-vector 

P-e A (P)=0, VA; e x (P) • e x ' (P) = S xx ' . (54) 

In ladder truncation, V2 = 0. From this point the solu- 
tion of Eq. (|49|) proceeds as described above for charged 
pseudoscalars. 

We are now in a position to fix the reference values of 
the mass-scale parameter and the current-quark masses. 



In the absence of electromagnetism we fix the values of 5 
and the sum (m u + ma) =■ 2m so as to obtain the exper- 
imental values of the ratio m w o /m p o and m p o . Moreover, 
the pseudoscalar variant of Eq. (|4"5)) produces degenerate 
bound-states associated with T 1 & T 2 - n^, two others 
associated with T A & T 5 - K^, another two with F e 

& T 1 - K°, K , etc. This is similarly true of the vec- 
tor equation: p ± , K*^, etc. We choose to fix the mass 
difference (m u — rrid) by requiring 

[m K ,o - m K , + ] f = 5.42 MeV, (55) 

which is the weighted average of the isospin-only differ- 
ences estimated in Ref. [HJ . The rainbow-ladder trunca- 
tion produces pure ss, cc and bb vector mesons and m SlC) & 
are set by identifying these states with the observed <fi-, 
J/ip- and T-mesons. This procedure yields the following 
values 

Q = 0.537 GeV, (56) 

m u md ms mc mb 

0.0140^ 0.0271^ 0.323 5 2.55 5 8.67 5 
7.51 MeV 14.6 MeV 173 MeV 1.37 GeV 4.65 GeV, 

(57) 

and they produce the meson masses in TablesU and [Til 
(NB. In solving BSEs to obtain masses the contribution 
from all orders in the current-quark mass splittings are 
incorporated.) We remark that the current-quark masses 
yield the following dynamical constituent-quark masses 
via Eq. |g7J (in GeV): 

M? M° Ml Ml M b D 

0.543 0.548 0.669 1.662 4.771 ' K ' 

For light quarks Mf - raj = Mq - m//4, where M<f 
is the chiral limit value [31(, and we note that in gen- 
eral M® — rn/ is a monotonically decreasing function of 
to/, bounded below by zero as to/ — > oo. This result em- 
phasises that the essentially dynamical component of chi- 
ral symmetry breaking decreases with increasing current- 
quark mass, as observed previously 0, H3|- 

The model we're employing is ultraviolet finite and 
the current-quark masses in Eq. (|57p cannot be directly 
compared with QCD's current-quark mass-scales. Never- 
theless, the values are quantitatively consistent with the 
pattern of flavour-dependence in the explicit chiral sym- 
metry breaking masses of QCD. It is notable that so far 
as isospin breaking is concerned, m u /md = 0.52, which 
is compatible with other contemporary estimates, e.g., 
Ref. 

Only £ and 9^ in Eqs. (|43f and ((44]) remain unknown. 
In order to fix these parameters we consider the neutral 
pseudoscalar mesons. In this instance Eq. (|4"9")l is modi- 
fied to the extent that Kl — > Kl + Ka', i.e., it reads 

r H (k;P)= [ [S{q+) 

Jq 

x T H (q + , q-)S(q-)] sr (K L + K A )%(q, k; P) . (59) 
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TABLE I: Vector meson masses calculated from the BSE de- 
fined by Eqs. ©, (g5J), (g3J), (O and (gHJ), using the parame- 
ter values in Eqs. (156[) and (157[1 . The experimental values are 
taken from Ref. [361 ]. The three parameters and current-quark 
masses were fitted as described in connection with Eqs. (|56[) 
and ()57[) . See Sec. lIV C~5l for further discussion of p° and w. 
|Expt. (GeV) Calc. (GeV) Th/Ex-1 (%) 





0.7755 


0.7704 


-0.66 


p ± 


0.7755 


0.7755 







0.7827 


0.7806 


-0.27 




0.8917 


0.8915 


-0.02 




0.8960 


0.8969 


0.10 




1.0195 


1.0195 





D*° 


2.0067 


1.8321 


-8.7 




2.0100 


1.8387 


-8.5 




2.1120 


1.9871 


-5.9 




3.0969 


3.0969 





B*± 




4.8543 




B*° 




4.8613 




Bf 




5.0191 




m ± 




6.2047 




T 


9.4603 


9.4603 






TABLE II: Pseudoscalar meson masses calculated from the 
BSE defined by Eqs. Q, (gU, flSJ, flU} and (gSJ, using the 
parameter values in Eqs. (|56l) and (|57|l . The experimental 
values are taken from Ref. [3g|. The three parameters and 
current-quark masses were fitted as described in connection 
with Eqs. ([56) and (1571) . 





Expt. (GeV) Calc. (GeV) Th/Ex-1 




0.13498 


0.13460 


-0.3 




0.13957 


0.13499 


-3.3 




0.49368 


0.41703 


-15.5 


K° 


0.49765 


0.42662 


-14.3 


V 


0.54751 


0.45499 


-16.9 


v' 


0.95778 


0.91960 


-4.0 


D° 


1.8645 


1.6195 


-13.1 


D ± 


1.8693 


1.6270 


-13.0 


Dt 


1.9682 


1.7938 


-8.9 




2.9804 


3.0171 


1.2 


B ± 


5.2790 


4.7747 


-9.6 


B° 


5.2794 


4.7819 


-9.4 


B° a 


5.3675 


4.9430 


-7.9 


Bt 


6.286 


6.1505 


-2.2 


Vb 


9.300 


9.4438 


1.5 



(%) 



As noted in connection with Eq. (|48|) . the gap equation 
is unmodified. It might be necessary here to empha- 
sise that this single equation describes all neutral pseu- 
doscalar bound-states; namely and for example, the only 
difference between the ir°, r\ and rj' is that they are ex- 
pressed by solutions of the BSE at different values of P 2 . 
From this perspective there is no mixing, as such, but 
one can pose the question: what is the quark content of 
each separated state. As we illustrate in Sec. lIV C~T1 the 
answers can be used to define mixing angles, all of which 
are in general different. 

The eigenvector for neutral pseudoscalars can be writ- 



ten 

d=0,3,8,... 

with the index selecting all Nt diagonal generators of 
U(Nf). Following the steps described above in connec- 
tion with Eqs. (j5"Tj) and (|52p . one can obtain the masses 
and Bethe-Salpeter amplitudes for the neutral bound- 
states. At this point we can fix £ and 0£ through a 
least-squares fit to the experimental values of and 
rrirj/mri'. This procedure yields 

£ = 0.076, 0^=0 (61) 

and the masses in Tablcllll 

C. Discussion of Results 

1. r/-ri' Mixing 

We'll begin with this topical issue and at first con- 
sider solving the BSE with an eigenvector of the form in 
Eq. (|60[) but with the sum running only over d = 0, 8; viz, 

r 80 (P) = 2F*g\P) + 2f°g°(P) , (62) 

with g' J (P) = 75 (ip\ + 7 • Pp 3 2 ) ■ This eigenvector splits off 
from that associated with T z and the ir° in the isospin 
symmetric limit. 

The case £ = provides a readily understood illustra- 
tion. In this case one obtains two bound-state solutions: 

mass (GeV) p\ pf, p ( { p 2 

m nn = 0.135 0.575 0.047 0.814 0.067 . (63) 

m Ss = 0.622 -0.786 -0.219 0.556 0.155 

Focusing on the eigenvectors, we rewrite Eq. (|62p in the 
form 

r 80 (0.133) = cos^.^ 8 ^ - sin 0i 2 J" ^ , 0i = -54.7°, 

(64) 

with g = 75(^1 + 7 ■ PP2) where p\ + p\ = 1. This is 
plainly a solution with ideal-mixing; namely, the light- 
est solution contains no s-quarks and is composed of an 
equal mixture of u- and d-quarks. The Bethe-Salpeter 
amplitude for the heaviest solution can be written 

r 80 (0.622) = sm9 2 2F 8 gl + cos9 2 2F°g 2 , 9 2 = -54.7°, 

which is a pure ss state. For £ = the dynamics de- 
couples from the flavour structure and hence the mixing 
angles for the two separated states are identical. 

With our preferred value of £ = 0.076, Eq. (fuTj) . we 
obtain 

mass (GeV) p\ p 2 p? p 2 

m n = 0.455 0.939 0.219 0.250 0.090 , (66) 

mw = 0.924 -0.260 -0.077 0.876 0.400 
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from which wc infer 



= -15.4°. 



V 



= -15.7° 



(67) 



Thus, while the Dirac structure of the rj and 77', described 
by <? 8 , 5°, is different, there is near equality between the 
mixing angle at each bound-state. For comparison, from 
a recent single mixing angle analysis one can extract [38[ 
9 = -13.3° ± 1.0°. The angles in Eq. (J6TJ) correspond to 
the flavour contents: 



\rf) ~ 0.55 (uu + dd)~ 0.63 ss . 
\t]') ~ 0A5{uu + dd) +0.78 ss . 

2. Chiral limit 



(68) 
(69) 



In the case of Nf — 3 massless quarks, Eqs. (0), (|4"5)l 
and (|59[) produce, without fine tuning, eight massless 
pseudoscalar mesons - the Goldstone modes - and one 
massive state. The massive state is solely associated with 
T° and 



m v > M =° 0.852 GeV, 



from which follows the model's value of 



(70) 



^L = (0.770 GeV) 2 



The chiral limit mass in Eq. ([70]) is 93% of the calculated 
value in Table ITT1 



3. -K^^q^q' Mixing 

With three flavours of quark, each with a different 
mass, all the neutral pseudoscalar mesons "mix"; i.e., 
there is no neutral pseudoscalar solution of Eq. (|59|) that 
is associated solely with a single generator of U(Nf). In 
this case the eigenvector assumes the form 

T H a (P) = 2 F 3 g 3 (P) + 2T s g s (P) + 2 T°g° (P) (72) 

and the BSE gives the solutions 

Pi 



mass 
(GeV) 
0.135 
0.455 
0.922 



Pi 



P2 



l>2 



P°l 



n° 
l>2 



0.002 0.009 0.001 
0.219 0.249 0.090 



0.996 0.081 0.023 
-0.026 -0.006 0.939 
-0.004 -0.001 -0.260 -0.077 0.876 0.400 

(73) 

From these Bethe-Salpeter amplitudes one infers the fol- 
lowing flavour contents: 



I''/) 



0.72 uu- 
0.53 uu 
0.44 uu ■ 



0.69 dd 
0.57 dd - 
0.45 dd- 



0.013 ss. 
0.63 ss , 
0.78 ss. 



(74) 
(75) 
(76) 



In the presence of a sensible amount of isospin breaking 
the 7T° is still predominantly characterised by J- 3 but 



there is a small admixture of ss. A glance at Eq. (|68|) 
shows that mixing with the tt° has a similarly modest 
impact on the flavour content of the 77 and rj ' . It's effect 
on their masses is far less. 



4- 7r°^7 Mixing 

There is merit in explicating the nature of the flavour- 
induced difference between the tt° and it masses. If wc 
ignore mixing with mesons containing other than u, d- 
quarks; viz., work solely within SU(Nf = 2), then the 
masses in Eq. (|57|) give ra T o — 771^+ — —0.04 MeV. On the 
other hand, it is apparent from Tables HI and ITT1 that the 
full calculation yields m^o — = —0.4 MeV, a factor 
of ten greater. When one considers only SU(Nf = 3); 
i.e., a so-called 3-8-mixing, then the 7r° mass is 0.1 MeV 
larger than in Table ITT1 m^a — ra^ = —0.3 MeV, and one 
obtains a mixing angle at the neutral pion mass shell of 



(77) 



For comparison, Ref. [39j infers a mixing angle of 
0.6° ± 0.3° from a X-matrix analysis pd — ► 3 He7r°. 
Plainly, mixing with the 77-meson is the dominant non- 



electromagnetic effect. Within this subspace, 



(71) 

larger than in Table |TT] and 



9^(m 2 ) = 1.3° 



is 5% 



(78) 



(Two mixing angles can be introduced to parametrise 
the complete problem of 7r°— 77— 77' mixing, e.g., (40j |. How- 
ever, that approach contains no information in addition 
to Eqs. tHI - 

It is noteworthy that 



W TO ?) ^ ^r,(m 2 o) 



inf 



' 7T7] 



= rl v e m] (ml ), (79) 
= 0.582 GeV -1 . (80) 



Our DSE framework describes mesons explicitly as 
bound-states of a dressed-quark and -antiquark. Hence, 
it is sensible to compare the result for this mixing an- 
gle and its momentum dependence with that, e.g., of 
Ref. [4l|]. They are commensurate: our mixing angle is 
< 20% smaller and the slope in Eq. |5D]) is < 20% larger. 
The slope in Eq. (|80|) is smaller than that which has been 
calculated in connection with p°-u mixin g b ut of the 
same order of magnitude; e.g., Refs. [42l. |43|. 44 1. 



5. Five flavours 

The tables present masses calculated with Nf = 5 
flavours of quark. The vector mesons are least compli- 
cated and hence we begin with them. It is a general fea- 
ture of the rainbow-ladder truncation that with no two 
current-quark masses equal, each neutral vector meson is 
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flavour-diagonal; namely, the kernel produces the follow- 
ing states, in order of increasing mass: uu, dd, ss, cc, bb. 
It is a good approximation for the heavier quarks and we 
therefore used this fact to fix the masses of the s-, c- and 
6-quarks. 

On the other hand, it is erroneous for the u, d-quark 
neutral vector mesons; viz., experimentally p° ^ uu and 
lu =/= dd, and that is why these states are written with 
quotation marks in Table [T] If one assumes m u = m<j, 
then there is no discernible problem. However, in the 
real case of m u ^ only an extended kernel can pro- 
duce the true flavour content for these states. Improve- 
ments along the lines pursued in Refs. [3(| HH, H2] do 
not ameliorate the situation because they preserve the 
flavour structure of the rainbow-ladder truncation. It is 
probable that inclusion in the kernel of diagrams which 
correspond to two- and three-pion intermediate states is 
required in order to provide the remedy, since the former 
contribute primarily in the uu — dd channel and the latter 
predominantly in the uu + dd channel. Solving the BSE 
without these channels but with forced (uu — dd) / y/2 and 
(uu + dd)/y/2 flavour contents, the p° and w are degen- 
erate with p ± . An analysis of diagrams corresponding 
to intermediate states containing pseudoscalar mesons is 
capable of lifting the p°-u) degeneracy, even for m u = 

Sa HE HE Sli - 

Regarding the other entries in Table [H where a deter- 
mination is possible the model is evidently accurate to 
< 10%. We therefore hold that it is reasonable to expect 
that the masses predicted for the as yet unobserved B*- 
mesons are similarly accurate. Indeed, they are probably 
an underestimate of the physical values by no more than 
this amount. 

One can read from Table Q] that for vector mesons 
composed solely of light-quarks the calculated flavour- 
dependent mass differences are well approximated by the 
differences in dynamical constituent-quark masses. The 
accuracy is better than 2%. We infer from this an es- 
timate of the non-electromagnetic part of the neutron- 
proton mass difference: 

[m„ - m p ]f = 5.3 MeV =: A UD . (81) 

This simple projection could be checked via the Fad- 
deev equation approach to nucleon structure [34J, which 
is kindred to that employed herein for mesons. For 
comparison, a numerical simulation of lattice regularised 
QCD has been used to infer a value for this difference 
of 2.26 ± 0.72 [49]; and the experimentally determined 
value, which includes electromagnetic effects, is 1.3 MeV. 

Additional context is provided by the observation that 
this simple reasoning also entails: Mj- — M^o = M^o — 
M s + = M s - - M 3 + = A UD . The experimental val- 
ues are, respectively: 3.3, 4.8, 6.5 MeV. We emphasise 
that pseudoscalar-meson self-energy diagrams contribute 
materially to a baryon's mass [5(| • When dealing with ef- 
fects of this small magnitude, differences in mass between 
the mesons that appear in such diagrams will contribute 
to these mass-differences. A complete, accurate calcula- 



tion as opposed to our estimates will naturally require 
precision. 

On the other hand, for vector mesons containing 
at least one heavy-quark, the difference in current- 
quark masses provides a better estimate of the non- 
electromagnetic mass difference. It is apparent from Ta- 
ble [IT] that the same is true of pseudoscalar mesons con- 
taining at least one heavy-quark. 

Owing to DCSB, the mass-squared of light-quark pseu- 
doscalar mesons rises linearly with current-quark mass. 
This explains our result that for the non-electromagnetic 
part 

[m K o-m K +]f = 9.6 MeV > [m K ,o-m K ,+] f = 5.4 MeV. 

It is noteworthy that for each system our calculated 
result for the H° - H+ mass-difference is larger in mag- 
nitude than the experimental difference. This is a neces- 
sary result and a useful check on this and other calcula- 
tions that omit electromagnetic effects. The inclusion of 
electromagnetism will act predominantly to increase the 
mass of the charged state and hence the mass difference 
will fall in magnitude. 

We have described above the quark-flavour content of 
neutral pseudoscalar mesons as calculated in subspaces 
of the full flavour group. This analysis can be repeated 
for U (Nf = 5) with the resulting flavour probability am- 
plitudes: 





uu 


dd 


ss 


cc 


bb 




0.72 


-0.69 


-0.014 


0.000 


0.0000 


\v) 


0.53 


0.57 


-0.63 


-0.022 


-0.0057 


W) 


0.44 


0.45 


0.78 


-0.060 


-0.0141 


\Vc) 


0.06 


0.06 


0.05 


0.995 


-0.0037 


M 


0.02 


0.02 


0.01 


0.005 


0.9996 



(83) 



As one might have anticipated, for states not much af- 
fected by Ka the probability of finding a particular type 
of "hidden flavour" drops as the mass of the quark flavour 
increases. This table indicates that the commonly used 
approximation of writing neutral light-flavour mesons in 
the form c\(uu + dd) + C2SS is accurate at the level of 
< 5%. 

6. Vector-pseudoscalar mass splitting 

As remarked in the Introduction, it is natural to ex- 
amine the manner by which the l~-0 _ mass splitting 
evolves with current-quark mass and mass difference. It 
is apparent from Fig. [2] that without exception the mass 
splitting, rfiH = mu* — rnn, decreases with increasing 
fhn = (1/2) (rriH* + inn)- The gross behaviour of the 
evolution is described well by a l/m#-dependence, which 
entails 

ra\, -m% ~ const. = 2p = 0.62 GeV 2 . (84) 

This outcome is consistent with observation. It is plain 
that since fhu is a measure of the dynamical constituent- 
quark mass, the global picture is not consistent with a 
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m H [GeV] 



FIG. 2: Vector-pseudoscalar meson mass splitting calculated 
from our results in Tables HI and HT1 namely, rfiH as a function 
of ffiH, where rhu = mw —ran and fhu = (l/2)(mfl« +mi). 
For the purpose of this figure, we compared <j> with the ss 
pseudoscalar described in connection with Eq. (I63|l . The solid 
curve is n/friH with /i = 0.309 GeV 2 . 

simple, single hyperfine interaction between constituent- 
like quarks. That is not too surprising owing to the Gold- 
stone boson nature of light pseudoscalar mesons. Our 
calculations show that such a picture only becomes rea- 
sonable for bound-states containing at least one c- or b- 
quark. 

V. SUMMARY 

In connection with pseudoscalar mesons, the axial- 
vector Ward-Takahashi identity is a powerful tool whose 
import should not be ignored. Following from this and 
resting upon the empirical observation that the rj' is not 
a Goldstone mode, we demonstrated exact chiral-limit 
relations that connect the dressed-quark propagator to 
the topological susceptibility. Furthermore, we extended 
the mass formulae derived in Refs. [1, [l]| to the case 
of electric-charge-neutral pseudoscalar mesons, for which 



flavour symmetry breaking entails non-ideal flavour con- 
tent. Our development confirms that in the chiral limit 
the rf mass is proportional to the matrix element which 
connects the rj' to the vacuum via the topological suscep- 
tibility. 

To illustrate the implications of the mass formulae we 
introduced an elementary dynamical model. This in- 
volved an Ansatz for that part of the Bethe-Salpeter ker- 
nel related to the non-Abelian anomaly which assumes 
the most general internally consistent form. It is a key 
and novel feature of our study that an anomaly contribu- 
tion is included within the Bethe-Salpeter kernel to yield 
meson masses. We thereby avoid the oft used expedient 
of enforcing anomaly constraints a posteriori at the level 
of a matrix of masses of unphysical-mesons. 

In addition to the current-quark masses our model 
involves only two parameters, one of which is a mass- 
scale. It was employed in a wide-ranging analysis of 
pseudoscalar- and vector-meson bound-states with an 
emphasis on the effects of SU(N f = 2) a nd SU(N f = 3) 
flavour symmetry breaking. Section FlV Cl details our find- 
ings, which are too numerous to recapitulate here. Suffice 
it to report that, despite its simplicity, the model proved 
elucidative and phenomenologically efficacious. Our re- 
sults, both their qualitative and quantitative aspects, 
should serve as a valuable guide in the future study of 
neutral pseudoscalar mesons using more realistic interac- 
tions. 
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